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1. Introduction 

Quark masses are parameters of QCD. They can be computed non-perturbatively on the lattice 
once an hadronic input is given. Sticking to the renormalization group invariant (RGI) definition of 
the quark mass, lattice results for the strange quark mass with 2 [|, |, g] and 2+1 [g, g] dynamical 
flavors differ among each other by more than 40 MeV. Such a spread can be ascribed to various 
systematics, mainly cutoff effects and the renormalization procedure, which in many cases relies 
on perturbation theory only. For the result we will present here the renormalization is performed in 
a fully non-perturbative way, while the lattice spacing is still varied in a small range only (roughly 
0.1 to 0.07 fm). 

We write the relation among the bare current (PCAC) quark mass m; and the RGI mass Mi as 



Mi = Z M (go)mi(go) 



(1.1) 



where i has to be interpreted as a flavor index. The factor Zm can be split into two parts, the first one 
connecting the bare mass to m,(/i), the renormalized one (in the Schrodinger functional scheme) at 
a given scale ji. 



mi(n) 



(1.2) 



Zp(g ,ajU) 

and the second one M/m(/i), which relates two renormalized masses and is therefore universal. In 
mass independent schemes (like the one we are going to describe) this second factor is also flavor 
independent [g]. Its computation is the main result we report about here. 

The renormalization group equations for the couplings are decoupled in mass independent 
schemes and read _ _ 

d £ or-s .A m i 



(1.3) 



with g(jlt) the renormalized coupling at the scale jU. The T and j3 functions can be expanded in 
perturbation theory 

-g 3 {b + b l g 2 + b 2 g A + ...} , z(g) g 2 {d + dig 2 + ...}, 



and the exact (i.e. valid beyond perturbation theory) solutions of eqs. ( |Q| ) can be written 



A = }i{bog 2 



Mi = mi(2b g 2 
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(1.4) 



(1.5) 



(1.6) 



where the RGI parameters A and M appear as integration constants. At the level of RGI parameters 
the connections among different schemes can be given in a simple and exact way. We thus regard 
them as the fundamental parameters of QCD. 



2. Running of the quark mass in the SF scheme 

In terms of Schrodinger functional (SF) correlation functions the bare PCAC mass reads 



1 (do + d )fA(xo)+c A ad*dofp(x ) 
2/p(*o) 



m(go, K) 



(2.1) 



x =T/2 
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where K is the hopping parameter. We use the plaquette gauge action and non-perturbatively 0(a) 
improved Wilson fermions [g]. The improvement coefficient ca is also set to its non-perturbative 
value from I 



]. For the definition of the SF with Wilson fermions we refer to QlOp and to Jl 1[ ] for any 
unexplained notation. Here we only recall that the correlation functions /a and fy involve matrix 
elements of the axial current and the pseudoscalar density operators, respectively. According to 
eq. ( |L2| ) the running quark mass at the scale jj, is obtained by multiplying the mass in eq. ( ^T| ) 
by the ratio ZA(go) /Zp(go,afj.). The factor Za is scale independent and can be fixed by Ward 
identities |[12|]. All the scale dependence in the running quark mass comes from Zp. In the SF it can 
be defined as 

Zp(g ,L/a) x ' Jl 



m - 



0, 0=0.5, T=2L 



(2.2) 



ML/2) ' 

where the boundary to boundary correlation function f\ takes care of the renormalization of the 
boundary quark fields and the constant c ensures Zp = 1 at tree level. Eq. (2/2) makes clear that 
in the finite volume scheme we are using, the normalization scale pL is identified with the (inverse) 
extent of the system, which, in turn, is uniquely fixed (up to cutoff effects) by the value of the 
running coupling g 2 {L). The non perturbative running of the quark mass is described by the step 
scaling function Gp(u) 



Op(k) 



m{n) = Z P (g ,2L/a) 
m{\i/2) a^o Zp(g Q ,L/a) 



(2.3) 



f(L)=u 



which can be viewed as an integrated form of the t function. Introducing also the step scaling 
function for the coupling o(u) = g 2 {L) (computed in [13]) and solving the following coupled 



recursions 



u Q = g 2 (L max ) =4.61 
o{u k+ i) = u k 



u k = g 2 (L k ), L k = 2 k L n 



w = 1 

wk = [nti Op(«i)] 



-l 



w k 



m(l/L n 



m{l/L k ) 



(2.4) 
(2.5) 



starting from the low energy scale 1/L max defined by uq = g 2 (L max ) = 4.61 to the higher scales 

m(l/L max ) 



m(l/L k ) 



For k > 6 contact can 



1 /L k , k = 0, 1, . . . , 8 (with Lq = L max ), one can obtain the ratio 
be made with perturbation theory, in such a way that the RGI parameters can be computed from 



eqs. (L5, L6) by using the 2-loop T function and 3-loop /3 function. 

In practice we calculated Op(u) at 6 different couplings in the range u = 3.33 . . .0.98, extrap- 
olating to the continuum limit lattice results from resolutions L/a = 6,8 and 12 (see |[T]] for details 
about the extrapolations). The change in scale by a factor 2 needed for a and Op is "easily" imple- 
mented in the SF scheme as it amounts to doubling L/a at fixed bare parameters (in other words 
we used lattices with L/a up to 24). We parameterize the results by the ansatz 



op(«) = 1 — ln(2)rfo« + p\u 2 + P2U 3 



(2.6) 



with free coefficients p\ and p2- Inserting it in the recursion eq. ( \2.5\ ), for the case k = 6, we obtain 

M 



1.297(16) at n = \/L mas . 



(2.7) 
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Using the result for L max A from [ 1 3 ] , we plot in figure [j] the running of the quark mass versus jU / A. 
In the plot we only include the errors from the coefficients w^; the uncertainties on A and M would 
simply amount to a shift of the axes. The figure shows that in this case perturbation theory works 
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Figure 1: The non-perturbative running of m. 2/3-loop refers to the 2-loop T-function and 3-loop j3- 
function, analogously 1/2-loop. 

surprisingly well down to rather small energies. This property is anyway specific of the quantity 
and the scheme considered here. As an example, the non-perturbative j8 function in the SF scheme 



shows larger deviations from perturbation theory at the lower end of energies reached here [13]. 

3. The strange quark mass 

In order to apply the described procedure to the computation of the strange quark mass we 
first need to calculate the factor Za(^o) /Zp(go,L mwi /a) for the bare couplings (i.e. lattice spacings) 
used in large volume simulations. As we are going to use results from the relevant set of j8 
values is 5.2, 5.29 and 5.4. Since this region is covered by the result for Za(^o) (and for Z A on , see 
figure ^) in [ 12], it remains to compute Zp(go,L max /a). At j8 = 5.2 the value is directly obtained by 
a simulation at L/a = 6 while at the other couplings we had to interpolate the results from different 
L/a (the appropriate one wouldn't have been an integer). The numbers are collected in table [j]. 



Clearly now Zp and Zm refer to the chosen discretization, whereas the ratio in eq. (2.7) is universal. 



p 


z P 


Zm 


5.20 


0.47876(47) 


1 .935(33)(24) 


5.29 


0.4936(34) 


1.979(25)(24) 


5.40 


0.4974(33) 


2.001(29X25) 



Table 1: Results for Zp and Zm for three coupling values. 

As we work with two degenerate flavors, instead of the strange quark mass what we actually 
compute is M re f associated with a kaon made from two degenerate quarks. We use data for ro(jc) 
and raps(fc) from [Q]. After extrapolating ro(xr) to the chiral limit (i.e to fc = K c ) we fit the product 
ro(K" c )mps(fc) as a function of K in order to determine fe re f defined such that 

(r (fQ> PS (K: ref )) 2 = (r m K ) 2 = 1.5736 , (3.1) 
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where we have used ro = 0.5 fm [|T^]. Finally we compute the PC AC masses at the bare parameters 
j3, K re f(j3) with j6 =5.2, 5.29 and 5.4 in volumes of approximately 1.5 fm. This reduces cutoff 
effects in the PC AC mass due to finite size [[If]]. The result for the RGI quantity is 

M ref = 72(3)(13)MeV, (3.2) 

the first error is statistical, the second is our estimate of cutoff effects obtained by comparing the 



result at the finest lattice spacing with that at the coarsest one. The number in eq. ( p.2| ) is consistent 



m l = ^(mi + +ml ) = (M + M S )B RG1 , (3.3) 



with the quenched estimate of the same quantity plj]. Assuming a weak dependence on Nf we 
re-interpret our result as an estimate in the 3-flavor theory (this assumption of course has to be 
checked in the future), and relate M re f to M s through the lowest order chiral perturbation theory 
formula 

1 

2 [ 

where M = 1 /2(M D + M d ), and Brgi is a constant of the chiral Lagrangian. For degenerate quarks 
of mass M re f eq. ( B.3[ ) reads m\ = 2M re fBRGi, which implies M re f = (M +M S ) /2. Using the relation 
M s /M = 24.4(1.5) we obtain 

M s ps 48/25M ref => M s = 138(5)(26) MeV . (3.4) 

Higher order contributions from chiral perturbation theory are expected to be around 10% and 
thus below the accuracy we have reached here. Eventually the result can be converted to the MS 



scheme at the scale 2 GeV by employing the 4-loop T and j3 functions [ ]18[ ] for Nf = 2. This yields 
m|p(2 GeV) = 97(22) MeV. 

We give our conclusions discussing the collection of results in figure ||[ By comparing differ- 
ent 0(a) improved regularizations we see that cutoff effects on the strange quark mass at a lattice 
spacing of 0.1 fm are around 20%. It is more difficult to assess the Nf dependence mainly because 
only perturbative renormalization has been used in the 3-flavor case. On the other hand the compar- 
ison of Nf = 2 determinations indicates that the use of perturbation theory for the renormalization 
constants underestimates the quark mass. 
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Figure 2: Summary of the strange quark mass data from lattice simulations. In the legend the discretizations 
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as suggested thereby. 
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